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The mechanism of /-electron delocalization is investigated within the multi-orbital Anderson 
lattice model by means of diagrammatic perturbation theory from the atomic limit. The derived 
, equations couple the intra-atomic transition energies, their spectral weights and population num- 

' bers of the many-electron states. Its self-consistent solution for praseodymium metal shows that 

04 ' the delocalization can be caused by external pressure via a resonant mixing of /- and conduction 

$_( , electrons in the vicinity of the the Fermi surface. It is also found that: 1. An increase of mixing 

leads to a decreeise of the physical values of the Hubbard interactions, U* , the reduction, however, is 
small. 2. The initial Hubbard U is split by renormalization into a set of different physical values of 
C/j*j. 3. The gain in cohesive energy together with the /-sum rule cause a transfer of spectral weight, 
' which is decisive for the delocalization of /-electrons. 4. The correlated Fermionic quasi particles 

have their bandwidth slightly reduced compared to the ones obtained by means of the Kohn-Sham 
equation. 

. 

. I. INTRODUCTION 

S: 

' Published as: Phys. Rev. B 62, 16370 (2000-11). 

^ The correct treatment of the local charge density built into the density functional method of ah initio band structure 
Q ■ calculations, particularly in the local (spin) density approximation (LDA), provides a remarkably good description 
Q ' of those materials where the energy structure of the ions, established by the intra-atomic Coulomb interactions, is 
strongly modified by the band formation. If this condition is not fulfilled, like in many materials containing d- or 
I ' /-electrons, the band structure method is insufficient. The parameter which determines if the material belongs to 
^ , this class of strongly correlated systems (SCES) is the ratio W/U, where W is the bandwidth and U is the parameter 
\^ • describing the Coulomb repulsion between d- or /-electrons. There are a number of different suggestions of how to 
CsJ ' treat the failure of the LDAjnethod. Also, at least three different definitions of the parameter U in use. The first one 
was introduced by HubbardEl as the Coulomb integral calculated between the wave functions of the corresponding d-, 
' or /-shells. Hubbard also discussed the importance of screening effects on U. This is mainly the Slater F° integral 
^—^ {Ui—F° is commonly used). The second definition, which comes from the formulation of the s-band Hubbard model, 
' is [/2 = [En-f i-En]-[En-Eri-i]7 whcrc E„ is the energy of a d(/)-atom with n electrons in the shell. An uncertainty 
. exists also within this definition. The recipe used in the so-called LDA+U methodo consists of calculation of the 
' energies E„_i and E„+i of an atom either within an impurity-type or super cell approach. The results are then 
^ ] depending on where, and from where, the electron is removed, the degree of relaxation allowed for the charge on the 
I ■ neighboring atoms and the number of atoms allowed to be involved in the relaxation. In the next step of the LDA+U 
' O method the found value U2 is used for correcting the potential experienced by d(/)-electrons only, in spite of the fact 
S that relaxation of charge on the neighboring sites has, been used when calculating U2- Another approach is to use 
^ U2 — S'^Elda/Sv?^ as discussed by Gunnarsson et al. a. 

• • , It is clear from a physical point of view that the intensity of the correlations characterized by the value of U should 
^ ' depend on the|-.type of ground state and the approximation used. Progress in this direction has been achieved by 
k>( ' Beiden et al. The definition used in their work, C/3 = /^^ + /^^ — 2/°^, involves both the Coulomb matrix elements 
5—1 linear muffin tin method using the spherical approximation (LMTO-ASA), , and the average population 

d ^ numbers of different shells, /"^ — [/°s?t-s + fppn^ + 2f^pnpns\/[np + ns]^. Both the matrix elements and the population 
numbers, are calculated self-consistently. The functional for these band structure calculations has been constructed 
on the basis of the Gutzwiller-type approximation, and been applied to the 3-d metals V, Mn, Ni, Cu and 4-d metals 
Nb and Tc. The Gutzwiller correlation correction influences the value of indirectly, via the variational procedure 
for finding the population numbers. 

If we look at the definition of C/2, the physical values of U* are determined by self-consistent values of the difference 
of the many electron energies E„. When we are discussing the ground state, or thermo-dynamical properties, all 
processes, fast and slow, renormalizing E„ have to be considered. The renormalization of A„^„_i=E„-E„_i by long- 
range Coulomb interactionQ leads to a dependence of these energies on the shape of the Fermi surface via the frequency- 
dependent dielectric function. Peaks in the electron density of states amplify the effect of mixing and, therefore, cause 
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a decrease of E„-E„_i, favoring a delocalization of (i(/)-electrons. In this paper we will consider renormalization of 
A„_„_i caused by kinematic interactions due to mixing. This will be considered in the present paper. The description 
of the high-energy experiments like photo-electron spectroscopies, sometimes requires a different approach (see Ref.Q). 

A very specific situation can arisejn metallic /-systems when one of the intra-atomic transitions happens to be 
in a close vicinity of the Fermi leveH. The scattering processes of the conduction electrons on this transition are 
resonantly amplified in this case. Therefore, one may expect a strong dependence of all renormalizing magnitudes on 
the mixing interaction if there is such a degeneracy. 

In this paper we show that the delocalization of /-electrons with pressure can be accounted for in a simple model. 
We also give quantitative numbers as regards the delocalization in Pr. Below, in order to display explicitly the 
mentioned effect of quasi-degeneracy we consider the periodical Anderson model for the case of praseodymium metal 
in the corresponding region of parameters. Here, by means of a diagrammatic perturbation theory from the atomic 
limit, we derive a system of self-consistent equations which determine the physical value of the Hubbard U, the spectral 
weights of the [{n — l),n] and [n, {n + 1)] electron transitions. The equations include correlation effects explicitly, 
without decoupling of the higher order correlation functions to the single-electron population numbers, and describe 
the dependence of the physical properties involved on the degree of delocalization of the states. As will be seen below, 
this procedure always leads to a splitting of the initial value for U to many different orbital-dependent U*. Here, 
we only inspect the orbitally polarized solution, the effect of crystal-field levels is not taken into account, since the 
scale of the crystal-field transitions is much smaller than the ones considered here, although it can be included in the 
formalism. 



II. THE MODEL AND ITS APPROXIMATION 



We start with the periodical multi orbital Anderson Hamiltonian, 

^ = E ^k,4^k, + E + ^ E U.,n,,n,, + E [^.(k)e*''^'4/,, + V;ik)e-^^^^ flc^^], (1) 

where n^^ = //^/i/x, i is the site index and fi is the /-orbital index, i.e. the cr, m; states. The first term describes the 
non-/conduction electrons, is the band structure energy, which in our case is taken from a self consistent band 
structure calculation. The two next terms describe the strongly correlated /-electrons, and the last term represents the 
mixing interaction between the /-states and the conduction band. Thus, we will study self-consistently the changes 
in the band structure caused by the mixing interaction in the background of the strongly correlated /-system. We 
define many-particle /-states, F, as 

|ro) =' |o) 



r^) = Im) = 4|o) 



and so on for all 14 /-orbitals. 

The creation and annihilation operators for the /-electrons are rewritten in terms of Hubbard A-operators, describ- 
ing transitions between many-particle states F and F', xf'^ |F}(F'|. The expansion is written as 



f^, - J2(^\f^^\^')xf^^ = ^(F|/,J/.F)Af (2) 
r,r' r 

We define the transitions F to /iF (F to F') as a, and the inverse, /iF to F, as a. Now let us reformulate the Hamiltonian 
and the Green functions (GF's) for the /-operators, in terms of Hubbard operators and then use a diagram technique 
for the X operatoraJ. The electronic Matsubara temperature /-Green functions should be expressed in terms of 
A-operator GF's 

= E {rl\U\^^r^){^'T,\fl,\^,) . {TXf^^^^^r)Xt'^^^^%r')). (3) 
ri,r2 

The X-operators satisfy the commutation relation 

[Xf■^ A^l = S., {S,rXP-^ + vSspX^-'^) . (4) 
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The fact that the (anti)commutator of two X operators is still an operator gives rise to "end" factors (averages of one 
X-operator, {X^^'^)) when the extended Wick's theorem is applied to the Green functions. In the diagram technique 
end factors are represented by open (closed) circles for the zero (dressed) spectral weights. The zero Green functions 
for X, and c-operators are 

(c^riA^Lo^))' = -*(rX«(r)xi(r'))[i = 5{a - b) . ^".^ W P>°, (5) 

(C^(k,zc.„))° -*(rck^(r)4(r'))° - , (6) 

lUJn - ek^ 

where iujn = (2n + 1)ttT and T is the temperature in energy units (k^^l). Pa denotes the end factors, Pa = Pv.v — 
{{X^'^\X^''^]) = {X^^^) + {X^' '^') = TVr + ^r' , which in lowest order are spectral weights, and Ar',r ^En^i-En^ 
the energy of the transition a entering the locatorgl3, Va{iuJn) — j_A ' ^'^^ X-operator Green functions. 
Averaging the expansion (||J^) for the anticommutator 

I r r' J r,r' 

gives a sum rule for the spectral weights in the X-operator Green functions. 

Using a diagram technique for the Hubbard operatorda we formulate a mean field theory (defined as the theory 
with no energy dependent corrections to the self energy, E, and to the spectral weights Pr,r')- The second order 

graphs are summed to give corrections to the population numbers, Pa = Pa + 5Pa^^ + 6Pa'^\ and the frequencies, 
Aa = Al + SAa- Within the Fermionic hierarchy, the first correction to the population numbers, the Ising-model-like, 
SPa^\ is given by the sum of ovals, and is shown in Fig.|^. 




FIG. 1. The correction SPi^^ to the population numbers. 



We denote the conduction Green functions by dashed lines, the locators for X-operator Green functions by full lines, 
bold full lines and double dashed lines mean the full Green functions, and the dots denote the terminal factors (local 
spectral weights). 

The-c,umulant contribution, SPa^^ (see Fig.|l|) can be calculated from a generating functional as described by Sandalov 
et al. M. 




(2) 

FIG. 2. The correction SPa to the population numbers. 

This procedure gives the expansion of an exponential function. However, this contribution, in the paramagnetic state, 
does not change the relative weights of the population numbers and can be omitted. The second contribution SPa , 
see Fig.^, is given by 

LJn ,/i,ai ,a2 
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where k are structure constants of the algebra and f3 ^ = T is the temperature and 



The Green function for the conduction electrons are dressed by the interaction as shown in Fig.^ 



(9) 



= = -->-- +--->- 



l-=>-- 



to give 



FIG. 3. The graph for the hill Green function for the conduction electrons. 



a 

The shift of the frequencies, SAa = Aa — A^, comes from the equation shown in Fig.^ 



(10) 




and is given by 



FIG. 4. Correction to the frequencies SAa. 



(11) 



These equations are coupled and accordingly have to be solved self consistently. Note that, so far, this formalism is 
not restricted to the case of praseodymium, it is valid for any /-electron system. The physics of /-systems is described 
by the positions of the levels and the corresponding spectral weights Pa- 

III. APPLICATION TO PRASEODYMIUM 

We now concentrate on the model for praseodymium metal. The existence of crystal-field effects says that Pr is in 
an orbitally polarized /^ state. Two localized /-electrons are placed into orbitals 1 and 2; this two-electron state is 
denoted as |12) (in indices as (12)). When mixing is absent n/ = 2 and N^i2) = 1- Switching the mixing on causes 
transitions to /^ and /^ states; we denote /^ as 1 or 2 depending on m;, /^ as |12z^), being one of the 12 unoccupied 
/-orbitals. The way that spectral weight transfers to higher transitions is shown in Fig.|^. 



1(12) 



(lv)(12v) 



(Iw )(12vv) 




(2v)(12v) 



1 2 body states 

2 3 body states 



■(12v)(12w) 



(2w')(12w) 3^4 body states 



FIG. 5. Possible transitions. The boxes indicate the states which have nonzero spectral weight in the start configuration. 
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The energies for these states are estimated from En = ne" + Un{n — l)/2 with e" providing the energy minimum for 
n = 2 and describing the attraction to the nucleus. This gives e° = —3U/2. The bare energies of the transitions are, 
from |12) to |1) or |2), denoted by indices as (12)1 and (12)2, Ai = E2 ~ Ei ^ A(i2)i = A(i2)2 = ~U/2, from |12i^) 
to |12) A2^E3-E2^ ^{i2u)(i2) = U/2 ^ A(i2^)(i^) = A(i2^)(2z.), and from \l2vv') to \\2v) A3 = 3C//2. The zero 
/-Green functions can be found by inversion of the expansion in Eqn. H, F^^ — {Tfif\)^ = F^'^ and Fy^^ = {T f^fl)^ 
becomes 



(0) 



pO 



AO 

^(12)2 



nO 

{I2v)(2v) 



AO 



- AO ' 
^(12)2 



using the expansion (Eqn, 



and for the 12 different v orbitals 



pO 

M12!^)(12) 



pO 

-^(12i/i/')(12i') 



AO 

(12>')(12) 



AO 

(12i/i/')(12iy') 



IUJ„ 



AO 

(12i')(12) 



Then, the dressed conduction Green functions ( [l0| ) can be written as, 

1 



Cp(k,ia;„) = 



|y^(k)|2F^(zu;„ 



(12) 



(13) 



(14) 



Due to the complex structure of the system of equations psaj make some further simplifications. We do not allow 
non-diagonal hopping (this case has been considered in ref.Eil) to occur. Non-diagonal terms affect the shifts only 
in fourth order of perturbation theory and are thus much smaller than the effects considered here. Physically the 
non-diagonal hopping would give a small width to the levels, not changing the physical picture to a large extent. 
Further we assume that the mixing parameter, V^(k), has a k dependence following ek- Then the summation over 
k in the formulas above can be replaced by a integration over the density of states (DOS) g^(e) ~ X^k"^!^ ~ ^^k)- 
Applying the formula (Eqn. 11) for the transition a=(12)l we obtain 



,5A(i2)l = l^2(k)pP(i2)l(ic^„)C2(k,ic^„) 



(15) 



and transforming summation over the Matsubara frequencies to a integration over energy, we obtain an expression in 
terms of the retarded Green functions 



'5A(i2)i = / dujf{uj-fi) / dxg2(x)\V2{x)\ 



I Im 



1 



UJ - 



id- A 



(12)1 



The correction to the spectral weight according to Eqn.( 0) becomes 



(12)1 



E 



dLuf{uj — jj) 



dxg2{x)\V2{x)\' ]Im 



Pi 



(12)2 



[uj + i6 ~ A(i2)2\ 



dujf{Lu — fi) / dxg,^{x)\Vi,{x)\ ) Im 



Pi 



(12iy)12 



_[uj + i6 - A(i2,.)(i2)]' 



C^{x,uj + id) 



C^{x,uj + i5) 



(16) 



(17) 



Note that the contributions from the lower and upper Hubbard sub-bands enter with different signs, i.e. these contri- 
butions favor a localization and a delocalization respectively. Due to symmetry of the corrections the (12)1 transition 
gives the same corrections as (12)2, and (12;^) (1;^) the same as {12v){2v). The corrections which we have taken into 
account are listed in tables || and ||. 



a 


5Aa 


(12)(1) 
{I2v){lv) 
(12;.)(12) 
(12;/!^') (12!/') 


(|V'2|"0(12)l(ia;„)C2>k 
{|V2pr'(12^)(l^)(»I.^n)C2>k 
{ Kpr'(12,.)(12) (i(^n)C„)k 

©(12i/'i/)(12i/')(it^ii)C''/)k 



TABLE I. Corrections to the frequencies. Note that only 11 of the 12 unoccupied 1/ orbitals (i.e. the i/' orbitals) contribute 



to the shift of the — transition. (Summation over k and integration over frequency is denoted as, ■ ■ ■ = (■ • •)k). 
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(12)1 
(12j/)(li.) 

(12!/)(12) 

(12iy:/')(12!^') 



5Pa 

(|Kr[©(12)2]'P(12)2Cl>k- 

— 12{|V,^| [I'(l2i.)(12)] -P(12i.)(12)C'^)k,i' 

(|V2p[©{12iy)(2i^)]^-P{12i/){2i^)C'2)k — 

— {|V1/| [I'(12i/)(12)] -P(12i/)(12)C'i.)k 

(|V^2|'[X'(12)l]'P(12)lC2)k + 
+ (|VlP[0(12)2]'P(12)2Cl)k 

ll(|K|^[P(12i/i/'){12i^')]^-P(12i'i'')(12i'')'^'')l< 



TABLE II. Corrections to the spectral weights, the factor 11 comes from 11 of the 12 unoccupied u orbitals (i.e. the v' 
orbitals). 
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The bare band structure of Pr metal consists of conduction bands, a local level, Ai, below the bottom of the 
conduction band, a level A2 slightly above the Fermi energy, ej?, and A3, which is much higher than ep- The input 
DOS is shown in Fig.|l^ (V=0). The transitions, Af^, from the crystal-field levels, | ± 1), are shghtly above the 
transition Ai of the singlet state and are empty in the ground state. We do not write them here explicitly, since 
the summation over the CP-levels gives the same for the lower spectral weight. Therefore, at T=0 K, this does not 
change the delocalization scenario. 

IV. ANALYSIS AND CONCLUSIONS 

The equations (|^)-(|l^) represent the general form of the mean- field equations for the multi-orbital Anderson model 
and some conclusions of geaeral natur&juan be drawn from them. Firsts as seen in the type of solution considered above, 
contrary to the GutzwilleiLJ approacholj, the description of excitations in terms of /-orbitals becomes meaningless since 
many-Fermion magnitudes arc involved and the proper way to describe electrons is in terms of collective atomic like 
excitations described by the Hubbard operators. The Fermi-like /-excitations are separated into different subgroups. 
The part described by the lower poles of the /-electron Green function (Fermionic transitions), is not important 
for thermodynamics. Inspection of the next corrections to these spectral weights shows that they depend on the 
spin dynamics. The upper-pole excitations form delocalized states, contributing to the cohesive energy, and the 
corresponding band structure depends mainly on the /-spectral weights in the upper pole. The latter being coupled by 
the sum rule to the lower pole, refiect the dynamics of localized spins and, therefore, includes in principle temperature 
effects into the band structure in a different way compared to usual approaches. Second^ the upper and lower 
frequencies (A2 and Ai) are moved due to renormalization by mixing interaction towards each other, giving a reduction 
of the effective U* . All other frequencies are also renormalized, thus leading to a splitting of the Hubbard U into 
several band structure dependent U* ^, = Ai"^^'"-* — A^?'" Third, when mixing is turned on, part of the electrons 

are placed into the mixed state (cj^^^f ''^ ), and the sum rule for the total spectral weight causes the occupation of the 
lower pole to decrease with the same amount, i.e. the system undergoes delocalization. Therefore, in the considered 
solution the value of the observable local moment turns out to be linked to the strength of the mixing interaction in 
the other energy region, near ep. Particularly, the localized moment observed in neutron experiments and the Curie 
constant in high temperature magnetic susceptibility should be proportional to the spectral weight of the lower pole. 
Note that this mechanism is difficult to extract from the one-electron GF's since the transfer of the spectral weight is 
hidden in the complex many-electron vertices, and cannot be obtained in a one-electron picture at all. This mechanism 
of delocalization does not work if the mixing interaction is self-consistently switched off by the correlations. 

The correction JP^^^ to the spectral weights contain the locator Va{iuj) in the second power, whereas the frequency 
shift SAa only in the first. Therefore, 6P'^^'> depends on the critical parameter Apr' — A* much stronger than the shift 
6Aa and it is the spectral weight transfer, not the decrease of the Hubbard U, which is decisive for the delocalization 
of /-electrons under applied pressure for the scenario considered. Increasing the mixing, under external pressure, 
causes spectral weight to transfer from the Ai- to the A2-energy region. 

Within this model, where hopping is not taken into account, the mixing of (spd)-states with the transition A2 
forms a pseudo-gap near A2, the DOS increases in the vicinity of and electrons fill in additional states, therefore 
the Fermi energy ep moves down. Thus, part of f- Hubbard's transitions influence the cohesive energy of the metal. 
Composing the /-electron's Green functions spectral weights from the X-operator Green functions, and taking small 
enough values of the mixing interaction, wc find that the same statement is valid also for f -electrons. Thus, the energy 
needed for this transfer is provided by the PV term in the free energy and by the contribution from /-electrons to 
the cohesive energy. 

At last, some conclusions specific to Pr can be added to what's given above. The deformation of the band structure 
also results in an increase of the effective mass, which should give an enhancemeat of the Pauli-part of the magnetic 
susceptibility besides the crystal-field mechanism described by White and FuldeEj. 

However, as described above, the value of the localized /-moment should be slightly reduced. Indeed, switching 

on mixing leads to the shifts: N^^^^ -> N^^^^ - SN^,^^ N^^^^ = ^ A^(i2,) ~ n,x E(cL^f ''^''"^) ^ «• This 
increases the spectral weight in the upper Hubbard band to 12ncx- Since the total spectral weight of G^j^g) equal 
to one, the spectral weight in a lower pole, at w ~ A(]^2)1j decreases by the same amount. 
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V. NUMERICAL CALCULATIONS FOR PRASEODYMIUM 



We made two numerical tests, one with a square shaped DOS using a mixing in the form Vk = Vi + V2(ek — co) 
and varying the constants Y\ and V2 and the position of the levels. The second type of numerical test were done for 
a DOS for the conduction band of the Pr metal. The DOS was obtained using a full-potential LMTO method with 



the two 4/-electrons in the core andJor the real dhcp structure (see FigjlO] V=0). In this case we used 14c = Vi only. 
According to Anderson's argumentli3 due to symmetry reasons the k-independent part of rtte single-electron mixing 
should vanish. There is, however, a local contribution to mixing from Coulomb interactionllj, 

^icVV/ = Vx{l\A!^!^\V)c^X-'-^ = Yc^X-^^^, (18) 

which is allowed to be non-zero. The non-local part we model by the term Y^i^v. — eo) providing Yi[e\^ — eg) = 0. 
Thus, we assume that the main contribution to the k-dispersion comes from the neighboring atoms which gives the 
k-dependence in both ek and Vk. We choose \J to be 10 eV. The system of equations was solved by means of a 
steepest descent method. The criterion for convergence used was that the sum of the squared difference of the input 
and output between two consecutive iterations, i.e. ^J^^Pa — SPa'^"''''^'^)^ + ^^i^^a — <5Ai'^°'^''')^ should be less than 
10^^. It is worth to mention that it is very difficult to achieve convergence near the transition from localized to 
delocalized and that we had to mix a very small portion of a previous iteration into the new solution. Further, a good 
start guess was important for the convergence. 

A. Constant DOS Case 

For the constant DOS we choose the bandwidth to be 8 eV. The population of the lower level for different cases of 
the mixing is shown in Fig.^. 



8 




FIG. 6. Population in the lower pole, when changing the parameters Vi and V2. 



From this figure we see that the population drops when the mixing is increased, but, there is also an interesting feature 
for V2 — 0.2 eV and above. When increasing Vi, there is actually a localization of the electrons. Since V2 comes from 
the bandstructure (Vi is the local mixing coming from the Coulomb interaction) we see that changes in structure or 
composition which reduce V2, can actually mean that we have a larger degree of delocalization. Hence, delocalization 
does not always come from increased mixing interaction, but decreased. The Hubbard {/-parameter was extracted 
and the result is shown in 
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O — o Vj=0.00eV 
A A Vj=0.20eV 



_^ Vj=0.40eV 



0,2 0,4 0,6 0,8 

(eV) 

FIG. 7. Effective Hubbard U, V , wlien clianging tlie parameters Vi and V2. 



Fig.^A and B, wliere Vi, and V2 is changed respectively. We see that U decrease with almost 10% when mixing is 
changed. This actually moves the upper transition closer to the Fermi-level, and therefore also affects the delocaliza- 
tion. There is, however, another important effect, the distance between the upper transition, A(i2)i, and the chemical 
potential 11. The effect of this is shown in Fig.|[ 
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FIG. 8. Population in lower pole, when changing the parameters A(i2)i — fi. 
In this figure, we see that one driving mechanism of delocalization is the parameter A(x2)i — /i- Further, we see that 



the behavior is almost independent for the different values of mixing. The best fit to the Pio^ {Pion, = Ni + N2) 
versus A(i2)i — /i is obtained from the function 



Plow ^ 



(A(i2)i - + 7' 



(19) 



Where a, (3 and 7 are constants with the values given in table III, 
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1.03465 


-0.0949886 




0.0551839 



TABLE III. Values of the constants determining Piow, in Eqn. y. 
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Using this equation we determine the value of A(i2)i — /i when complete delocalization occurs to be 0.0366 eV. i.e. 
when the transition is very close to the chemical potential. 



B. LMTO-DOS Case 



When we used the DOS from a band-structure calculation we obtained the following result. For this case we 
put V2 — 0. The DOS after the self-consistent cycle is shown in Figjl^. As seen, upon an increase of mixing, the 
bandwidth is narrowed and a gap where the upper level is situated is opened up. The states, pushed out of the gap 
region, enhance the DOS at the Fermi level. The gap that develops and the delocalization is further illustrated in 
Fig.|l^. The peak above the Fermi level will not be seen in experiments since these transitions become delocalized due 
to mixing interactions. 
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FIG. 9. Changes in A'^(i2), U and fi for different choices of parameters. (V2 = 0). 
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V-0 



FIG. 10. The delocalization and changes in DOS for different values of the mixing parameter. {V2 = 0, and V = Vi in the 
figure) . 



In Fig.^ we present the result of the self-consistent calculations. As seen from equations (^6|) and (0), the energy 
difference between the upper level and the chemical potential is one of the most crucial parameters of the model. 
Fig.^ shows the dependence of the population number N12 on this parameter for two values of mixing. The Hubbard 
parameter chosen is 10 eV and the mixing, Vk = 14, is 1 and 2 eV respectively; the bandwidth is 9 eV. The next figure, 
Fig.p^, shows how U, defined as A(i2y)i2 — A(i2)i, change as a function of the same parameters. The graph in Fig.||d 
shows the change of the chemical potential /i and the fourth graph, Fig.^, shows the change in U and population 
number iVi2 when only mixing is changed. As seen, in Figj^c, dependences on the value of the mixing parameter 
are smooth and the only instability is caused by the level slightly above the Fermi energy. The value of the critical 
parameter A|'-|^2)(i2,y) ~ M can be changed by chemical doping, which will shift the chemical potential and therefore 
affect the transition pressure. If we look at the form of the DOS, generated by pure LDA calculations (see Fig.|l^ for 
V=0), it does not contain a quasi-gap in the vicinity of ej, since within band structure calculations all /-states are 
taken in core and the transition A2 does not exist. It is, however, worth to note that experimentally photo-electron 
spectroscopy shows peaks above the Fermi level for all rare earths. These peaks are interpreted as /-states. If we 



treat the /-states as fu 
described (see e. g. ref. 



delocalized we have peaks above the Fermi-level but the equilibrium properties are not well 



3). 



VI. CONCLUSION 



We have studied the scenario of /-electron delocalization and dependence of Hubbard U on the value of the mixing 
interaction which is assumed to increase proportionally to the external pressure. The scenario considered emphasizes 
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that the dominant role is played by the /-sum rule, connecting excitations in different energy regions and resonant 
mixing of (sprf)-clcctrons and Hubbard's /-excitations. The mechanism can only exist due to strong correlations. 
Thus, in our picture, delocalization arises due to the following reasons; First, the strong correlations provide the 
separation of the /-shell into two manifolds in different energy regions. Second, under the application of pressure the 
mixing increases due to an increase in the overlap. Third, the shift of bands, and the renormalization of the levels, 
alter the crucial parameter A(i2)i — /U. 
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